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Abstract 

We show how to reduce the non abelian Born-Infeld action describ- 
ing the interaction of two D-particles to the sum of elliptic integrals 
depending on simple kinematic invariants. This representation gives 
explicitly all a' corrections to D-particle dynamics. The a' corrections 
induce a stabilization of the classical trajectories such as the "eikonal" 
which are unstable within the Yang-Mills approximation. 



1 Introduction 



Since Dirichlet-branes (D-branes) found their true place in string theory [|l|], 
they have been challenging our basic intuitions about space-time. One of 
the most intriguing features is the fact that the coordinates describing the 
relative positions of D-branes naturally appear to be matrix valued 0, a fact 
that is at the heart of the matrix theory conjecture pi (for a review see 

In the matrix theory limit J5], ||, the form of the action relevant for the 
dynamics is simply the dimensionally reduced maximally supersymmetric 
Yang-Mills theory (7|- but this is not so in other regimes. A specific 
example is ordinary weakly coupled type II (A or B) string theory at fixed 
value of a', where the relevant effective action is the disk generating func- 
tional, already a highly non local object that is not known in closed form. 
If one restricts oneself further to considering the limit where the accelera- 
tions and higher time derivatives are small (keeping however the possibility 
of having relativistic velocities), the relevant action is the non abelian gen- 
eralization of the Born-Infeld action (NBI) [14[] to be described below[|. The 
full supersymmetric extension of such action is still out of reach, but for our 



purposes the bosonic action of |14) will suffice. 

One of the hurdles in using the usual form of the NBI lagrangian is the 
presence of the symmetrized trace operator in front of the usual square root of 
the determinant. This is particularly disconcerting in the case of a D-particle, 
where one would like to use this lagrangian to investigate corrections to the 



scattering processes of [|13[] or to the canonical quantization of [|TT], p^| . 

In this paper we solve this problem by showing how, in the case of two D- 
particles, the NBI lagrangian can be recast in the form of an ordinary function 
of a small number of kinematic invariants. The functional dependence of the 
lagrangian on these invariants in through elliptic integrals. 

As a first application of our result we investigate the stability of the 
eikonal trajectory, which is known to be unstable in the Yang-Mills system 
which is chaotic P^fl . We show how the a' corrections change the problem 
from the usual one of particles in flat space subjected to a quartic potential 
with flat directions to that of particles in conformally flat space with a similar 
potential. The curvature computed from the conformal factor is positive in 
the region near the flat directions thus contributing to the stabilization. 



1 For earlier references, see ifiE 
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Two further issues that we leave for future work are the possibility of 
existence of trajectories other than the eikonal but with the same asymptotic 
behavior as \t\ — > oo. It would be interesting to find such trajectories where 
two D-particles come close, exchange their identity through a rotation in 
"non commutative space" and then separate again. Such trajectories do not 



exist in the truncated Yang-Mills theory |25[ but may be present for the 
full NBI action. Another interesting issue would be to study the effects that 
the a' corrections have on the Born-Oppenheimer approximation, perhaps 
suggesting a change of variables that makes it applicable for small distances. 
The paper is organized as follows. In the next section we state Tseytlin's 



result [14| on the form of the NBI action that is the starting point of our 



investigation. We also summarize the basic elementary features of the sym- 
metric trace necessary to give a precise meaning to Tseytlin's action. In 
section three we show how to rewrite the NBI action as a simple function 
of a few invariants in the case of two D-particles. We begin by stating the 
result and all the assumptions that go into it and end with the actual proof. 
The last section contains a first look at the dynamics that can be obtained 
from such an action (e.g. the stabilization of the eikonal) and some possible 
future directions to be explored. 



2 The Non Abelian Born Infeld Action 

Already at tree level in the string coupling g s (disk diagram) and in flat 
space-time g^ u = 77^ = diag(— 1, 1, • • ■ , 1), the effective action for iV Dp- 
branes is a highly complicated object that is not known in closed form. This 
is of course due to the fact that, even at tree level, the n-point function 
on the disk receives contributions from the massive string states that, when 
integrated out, yield a non local functional of the massless modes. 

If one makes the further approximation of neglecting higher derivatives of 
the field strength F and, by consistency, terms involving the commutator of 
two field strengths as [F, F] = [D, D]F, one can cast the bosonic part of the 
remaining action in the non abelian generalization of the Born Infeld action 
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(NBI)B 

S = ~ J dxP +1 Str ^-det (^ + 27ra'F M ^, (1) 

where v = 0, • • ■ , 9 but the field strength F^ u only depends on the first 
p + 1 coordinates (coordinates on the Dp-brane). The determinant det is 
taken only on the Lorentz indices /i and v and the symmetric trace Str will 
be discussed in more details below. Note that the determinant of a matrix 
with non commuting elements can be defined in many inequivalent ways; the 
symmetric trace Str picks out the definition of physical interest. 



2.1 The Symmetric Trace 

The symmetric trace of a set of n matrices is defined as the sum of the 
ordinary traces over all possible permutations n G S n with the appropriate 
weight 1/n! 

Str(Ai,--.,i4 n ) = - £ Tt^d-.-A^)). (2) 

The trivial technical point that must be understood about the symmetric 
trace Str is that, contrary to the more familiar trace operator, it does not 
allow one to perform the matrix algebra inside it. To give an explicit example, 
let A, B and C be three arbitrary matrices and let D = AB. Then, while 
obviously Tr (ABC) = Tt(DC), one has from (^) 

Str(A, B,C) = ~ Tr(ABC + ACB) ^ Str(Z>, C) = Tr(ABC). (3) 

That's why, to avoid any confusion it is better to write Str(A, B, C) instead 
of Str (ABC). For an arbitrary function f(A±, • • • , A n ) of n matrices Aj, the 
symmetric trace (with respect to those matrices) is defined by first expand- 
ing / as a Taylor series and then performing the symmetric trace on each 
monomial: 

Str(/(A 1; ■ • • , A n )) = £ 1 /fr."*J(0, • ■ ■ , 0) Str(4 fel) , ■ • ■ , A^). 

(4) 

2 In all other equations we shall drop the dependence on g s and 2ira' 
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It is important to understand that the superscripts (ki) inside the Str oper- 
ator simply mean that the matrix A4 is repeated ki times in the list and do 
not represent a multiplication. 

In the case of ([!]) by Str we always refer to the symmetrized trace with 
respect to the form F. 

3 NBI action of D-particles 

In this section we derive the explicit expression (i.e. we perform the sym- 
metric trace) for the NBI action of two D-particles in a particular gauge that 
although simplifies some of the computation, does not make us loose gener- 
ality. Since the computations involved are rather laborious, we will begin by 
explaining the choice of gauge and stating the result; the remainder of the 
section deals only with the derivation of our expression for the NBI action. 

3.1 Choice of gauge and statement of the result 

Let us consider two D-particles moving in the 1, 2 plane and let us set 
A = A, A k = X k (k — 1,2 only) so that 

F ok =j t X k -i[A,X k ]=DX k and F 12 = -i[X u X 2 ] (5) 

where all fields are SU(2) valued. While keeping track of the non abelian 
nature of the coordinates, one can still partially fix the gauge to beQ 

X k = x\a l + x\a 2 , and A = acr 3 . (6) 

This gauge can be reached as follows: at each time t one chooses the 
gauge transformation g(t) that sets the commutator [Xi, X2] proportional to 
a 3 . The remaining U(l) gauge freedom of rotating along the a 1 , a 2 plane 
is compensated by the gauge field A = aa 3 . If the commutator is zero, the 
transformation above is ill defined, but in this case X\ and X 2 are parallel in 
the SU (2) Lie algebra and can both be made independent on a 3 by a gauge 
rotation. 

3 Upper indexes always refer to Lie algebra components, lower indexes always refer to 
Lorenz components. 
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The main statement is that the NBI action in this gauge can be written 
in terms of the following three quantities, invariant in the unbroken gaugejj 

4(^1^2 — x 2 x 1 ) 

1 



p 

V 



X 



1\2 



+ (±l) 2 + (xiy + a 



2\2 



™2\2 



X 



1\2 



2 + (xlY 



{x\Y - {xlY 



+ a(x\x\ + x\ 



2 X 2 



2\2 



;,2 



) 2 ) 



(7) 



as 



S{p,T),v) 



-1/(1 + P) 



1 



:(F((j),m) - F(n/2,m] 



l-r) V z/ + P 
+ v ^Tp(£(^, m) - £(tt/2, to)) (8) 

where F and are the elliptic integrals of first and second kind respectively 
whose arguments are: 

, . /l — v 77 + p 

<p = arcsm * / and m = . 9 

V V l-r] v + p V ; 

Since there are different conventions in the literature, let us also recall the 
definition of the elliptic integrals we are using: 

d6 



and EU,m)= [ Jl-m sin 2 9 d9. (10) 
J o 



F W>> m ) / R ^TE 

Jo v 1 — msir 6 

F(tt/2, m) and E(7r/2,m) are known as complete elliptic integrals and are 
sometimes indicated by K(m) and E(m). 

For completeness, let us write the action of R, the generator of rotation 
in the 1, 2 plane, and G, the generator of the left over U(l) gauge invariance. 



R 



x\ 



, 2 

.A.- 

, 2 

.A.- '-) 



cos — sin 9x\ 
sm.6x\ + cos 

cos 9x\ — sin 9x\ 
sin #x 2 + cos 9x\ 



(11) 



For convenience we indicate the covariant time derivative of the components by a dot: 
DX k = x\a x +± 2 k <j 2 . 
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G : < 



cos a(t)x\ — sma{t)x\ 
sma(t)x\ + cos a{t)x\ 

cosa^)^ — sma{t)x 2 
sin a{t)x\ + cos a{t)x 2 



a — 



Ida 
2 dt 



It is an easy matter to check that (0) are invariant under R and G. 



(12) 



3.2 Derivation of the form of the action 

The remainder if this section contains the derivation of (|8|). 

We begin by evaluating the determinant in (|l|) without performing the 
matrix algebra:^ 

/ -1 



det(?7^ + F^u) = det 



DXx 



DX 2 
-i[Xi,X 2i 
1 



-DX X 1 
\-DX 2 i[X 1 ,X 2 ] 

= -1 + (DXi, -DX 2 , -i[Xx,X 2 \) + (-DX h DX 2 ,i[X h X 2 \) 

-(1,DX 2 ,-DX 2 ) - (1, DXi, —DXx) ~ (-hi[X 1: X 2 ],-i[X u X 2 }) 

= -1 + (DX k ,DX k ) + ([Xx,X 2 ], [X h X 2 ]) (13) 

where the last simplification is allowed because all these terms will be used 
inside the symmetric trace. By expanding the square root as a double power 
series one gets 

1 



S 



Str J - det (rj^ + F f 



2 

--Stiv I 



/IV J 



{DX kl DX k )-{[X ll X 2 ],[X 1 ,X 2 ]) 



m.q 



1 (2m + 2g-2)! 
2 2m + 2 ? m\q\{m + q - 1)! 



E 



Str( J DXf m) ,[X 1 ,X 2 ]^) 



(14) 



where, for m = q = we define the coefficient of the symmetric trace to be 
equal to its "analytically continued" value —1/2. 

5 The meaning of symbols like (DX\, —DX2, —i[Xi, X2]) is simply that of an unordered 
list of matrices, just as in (Q). 
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Eq. ( |H| ) can be simplified by noticing that the commutator [Xi, X2] anti- 
commutes with the covariant derivatives and has square proportional to the 
identity matrix. This allows one to write 

S«(DXr, [ Xl ,X4<*>) = I |^(^±*(-l)VStr( M r') (15) 

in terms of the invariant p defined in (^). Note that we are still using Str. 

The reader can check the validity of by chosing the term in the 
sum containing only DX\, expanding the symmetric trace into (2m + 2q)\ 
ordinary traces according to (§) and grouping the 2m\2q\ identical terms 
coming from the permutations of DX\ with themselves and [Xi,^] with 
themselves. In each of the remaining (2m + 2g)!/2m!2g! terms one then 
brings all the commutators [Jfi,^] to the left keeping track of the signs. 
There is an excess of (m + q)\/m\q\ terms with positive sign, yielding 

Eq. (|T5p follows by replacing (DX^ by {DX^ and Tr by Str on the 
r.h.s. of (H). 

With this simplification: 

n^=o 2 q m! <? ! 2 2g + 2m- V y 1 ' y ' 

To proceed further, define three (non invariant) intermediate quantities 

a = {x\f + {x\f 

c = {x\f + (xi) 2 (18) 
and carefully expand the symmetric trace using (0) 

Str( J DXf m) ) = Str (((xW + x\a 2 )^ + (d^ 1 + A 2 a 2 )( 2 )) M ) 
= Str (L{<j\ a 1 ) + b(a\ a 2 ) + c(a 2 , a 2 )) M ) 
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i+j+k=m ^'3' 

mPaW^_ / l + (-l)A (2i + j)\(2k + j)\ 
i+ fy =m 2m! i\j\k\ [ 2 ) (i+j/2)\(k + j/2)V [) 

The reason for the combinatorics in the last line of (^) is exactly the same 
as the one given for (15); also note that the symmetric trace vanishes for 
j odd explaining the presence of the projection operator (1 + (— l) J, )/2 and 
term of the type j/2 in the factorial. In fact, by using the identity 

oo oo 

E E a m,i,i,k = E a m+i+j,i,j,m (20) 
m=0 i+j+k=m i,j,m=0 

and replacing j by 2j, we have the following multiple power series for the 
NBI action: 

s _ ~ (-1)" 2q\ (2i + 2j)!(2m + 2j)! pV6^c m 

~ q ~t j=Q 2 2m + 2i +V+ 2 4 i\2j\m\(i + j)\(m + j)\ 2m + 2i + Aj + 2q-l' 

(21) 

To simplify things slightly, let us introduce four new variables x, y, z, w and 
define 

~ 2q\ (2i + 2j)\(2m + 2j)\ x 2m y 2i z ij w 2q 
U{x,y J z J w) = 2^ 



, J=0 q\ 2 i\2j\m\(i + j)\(m + j)\ 2m + 2i + 4j + 2q - 1 

(22) 

Trivially, if one knows U, the NBI action is 

S=U{^/2,^/2,Vb/2,i y /p/2). (23) 

To obtain a closed expression for U notice that there is a very similar 
expression that can be easily summed: 

7( , ^ 2g! {2i + 2j)\{2m + 2j)\ 2m 2l 4j 2q 

Z(x,y,z,w)= 2^ — x y z^w q 

q,m£j=o iml (z + j) ! {m + j) ! 

' (24) 
1 -4w 2 )((l -Ax 2 )(l -Ay 2 ) -4z 4 
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The functions U and Z are related by the PDE: 

(xd x + yd y + zd z + wd w — ljU = Z 



(25) 



and we are left with the problem of finding the solution to (|25|) compatible 
with the conditions given by the power series fl22]) . Let us first find the most 
general solution to (EST) by making the change of variable 



p = (xyzw) 1/A , r = ( — 



^ [v. 



1/4 



s I — 



(26) 



that reduces the PDE to an ODE 

(pd p -l)U = Z (27) 
whose most general solution in terms of an arbitrary function g(r,q,s) is 



U = p (g(r,q,s) + J ^ dp\ . 



(28) 



Performing the elliptic integral we obtaii 



U 



1 - 2(x 2 + y 2 + J(x 2 - y 2 ) 2 + z 4 ) 

P9{r,q,s)- ) =H1 

\ 1 - 2(x 2 + y 2 - ^J(x 2 -y 2 ) 2 + z 4 



F(</>, m) 



4w 4 ) 



2{x 2 + y 2 + J(x 2 -y 2 ) 2 + z 



2(x 2 + y 2 + ^J(x 2 - y 2 ) 2 + z 4 ) - 4w 2 



+ y 2(x 2 + y 2 + ^(x 2 - y 2 ) 2 + z 4 ) - 4w 2 £(0,m) 



(29) 



with 



arcsm 



1 - 2{x 2 + y 2 + ^(a; 2 - y 2 ) 2 + ? 



^ 1 - 2[x 2 + y 2 - 



y 2 ) 2 + z L 



m 



2(x 2 + y 2 - ^(x 2 - y 2 ) 2 + z 4 ) - 4w 2 
2(x 2 + y 2 + ^(x 2 - y 2 ) 2 + z 4 ) - Aw 2 



(30) 



3 the definition of the elliptic integrals is given in (fLCM) 



The function g(r, q, s) can be fixed by imposing that the solution should 
not have any term linear in x, y, z, w in its power series around the origin. 
Hence, by computing U along the line x = t/r 4 , y = t/q 4 , z = t/s 4 , w = t 
and imposing the vanishing of the 0(t) terms one obtains 

2(x 2 + y 2 + J(x 2 - y 2 ) 2 + z 4 ) 
Pg{r,q,s) = =F(7r/2,m) 

J2[x 2 + y 2 + ^/(x 2 - y 2 ) 2 + z 4 ) - Aw 2 

-^2{x 2 + y 2 + ^J(x 2 - y 2 ) 2 + z 4 ) - Aw 2 E(tt/2, m). (31) 

Tracing back the various changes of variables to (0) (c.f.r. also (|18|) 
and fl2"3"D we obtain our final answer (^). Note that the three non invariant 
quantities a, h an c have combined into two invariants v and 77 (p was invariant 
from the beginning) so that the final expression @ is depending only on the 
three variables (0) and not four. 



4 Some simple properties of the action 

To get a first feeling for the quantities p, i] and u, note that, in the commuting 
limit p = 7] = and v = v 2 , v being the relative velocity of the two D- 
particles. Substituting these values into (||), one finds the expected result 

5 = -v / T^= -v / T^2. (32) 

Another trivial limit is that where all time dependence is dropped, i.e. rj = 
v = 0, p = — Tr([Xi, X2} 2 ) where one finds the potential 

S = -yjl+~p = -~ Tr y/l - [X u X 2 } 2 . (33) 

With the proper factors of the string slope reinstated, to leading order in a' 
(|33|) is simply the Yang-Mills potential and p measures the distance from the 
bottom of the valley. 

A less trivial limit comes by asking what is the form of the kinetic term 
exact to all orders in a'. Here one already encounters an expression that 
could not have been guessed simply from (p. Although it is possible to 
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obtain this from the full answer (|8|), we find it easier to go back to the partial 
sum ([H]) and specialize to the m = 1 case. The coefficient of 

(2) \ — rp^f nv 2\ _ off -1\2 | /-1\2 | /-2\2 | /-2\2 



Str(DXn = Tr(DXf) = 2[{x\) 2 + {x\f + {x\f + {x\f) (34) 
is 

f 1 2„! (-1)9 arcsinh(yp) 

Thus, for small velocities but to all orders in a', the lagrangian governing 
the dynamics of two D-particles in this gauge is given by 



arcsinh ( 2{x\x\ — x\x\) , 
£ = TtV^ T^T-^fe) 2 +(i2) 2 + (^) 2 + (^) 2 ) 



-yjl +A(x\x 2 2 -x\x\) 2 . (36) 

Eq. (|36| ) describes the motion in a conformally flat space. 

It is possible to use the same change of variables introduced in |]24| to 
eliminate two of the four variables in (^) by virtue of the symmetries R and 
G. Setting 

x\ = f cos 9 cos (p — g sin 9 sin <fi 
x\ = — /cos^sin0 — g sin 9 cos (f> 
x\ = f sin 9 cos <p + g cos 9 sin <f> 

x\ = — f sin9sm cf> + gcos9 cos (p (37) 

and eliminating the cyclic coordinates 9 and <fi through their equations of 
motionF] 

(f + g 2 )4, + 2fg8 = 

Vrf+W+M = <— ^ (38) 
one obtains the Routhian (still denoted by C) 

C=\K (f 2 + g 2 ) - V, (39) 

7 l is a constant of motion corresponding to the orbital angular momentum. 
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where 



arcsinh(2/ff) I 2 f 2 + g 2 r 

K= — and V=— IF —^ + ^l+4f 2 g 2 . (40) 

In these coordinates the eikonal takes the form 



= 0, = 0, / = VvH 2 + b 2 , 6 = arctan(wt/6), (41) 

where v and b are the relative velocity and impact parameter respectively. 

It is interesting to compute the curvature associated to the conformal 
factor K. This is partly because a previous analysis |24| has shown that the 



eikonal trajectories are in fact unstable in the Yang-Mills case, the system 
being chaotic. A positive curvature would have a stabilizing effect on the 
classical trajectories. The Gauss curvature coming from the conformal factor 
K turns out to be 

R= . 2( f + g ^J^l, (42) 

and it is indeed positive near the eikonal, thus having a balancing effect on 
the small fluctuations around it. 

To study the balancing effect we analyze the system fl39|) numerically 
using Poincare sections. The Poincare section is a surface in phase space. 
Our phase space is four dimensional but energy conservation restricts the 
motion to be on a three dimensional subspace. One can then impose any 
constraint and fix a two dimensional surface named Poincare section. We use 
the constraint g = 0. A point is plotted each time the trajectory goes through 
the surface. A chaotic distribution of points says that the system under 
consideration is chaotic. If there are solid or dotted lines on the Poinacre 
section the system is regular. 

We perform about 1000 experiments for 10 different values of the energy. 
We plot the typical Poincare sections to show an appearance of the regular 
motion. Figures [1] and |] show the Poincare sections for trajectories started 
at point / = 5,g = and / = 5, g = 0.01, respectively. Comparing these 
figures we see that the eikonal trajectory is stable. Figures |3| and [| repre- 
sent stochastic behavior of trajectories with the initial coordinates / = 5, 
g = 0.1 and / = 5, g = 0.2, respectively. All trajectories above have zero 
initial velocities. Figures [| and |^ show the motion of the particle with initial 
conditions / = 10, g = 0.01, / = -1 and / = 20, g = 0.001, / = -1. One 
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can see that the particle reaches a point of minimal /, "scatters" and turns 
back without stochastization. 

It would be interesting to pursue the study of (Q) or some of its simplified 
limits like (|3~9p further, in particular addressing some of the issues raised in 
the introduction such as the detail study of trajectories that resemble the 
eikonal only for \t\ — > oo or ways of improving on the Born-Oppenheimer 
approximation near p = 0, interesting for the study of Matrix black holes 



It - 2 
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Figure 1: Poincare section for trajectory with initial data / = 5,g = 0. 
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Figure 2: Poincare section for trajectory with initial data f — 5,g — 0.01. 
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Figure 3: Poincare section for trajectory with initial data / = 5,g — 0.1. 
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Figure 4: Poincare section for trajectory with initial data / = 5,g — 0.2. 
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Figure 5: Poincare section for trajectory with initial data / = 10, g = 0.01 
and / = —1. 
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Figure 6: Poincare section for trajectory with initial data / = 20, g = 0.001 
and f — —1. 
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